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ABSTRACT 

We use high-resolution cosmological hydrodynamic simulations to study the angular momentum 
acquisition of gaseous halos around Milky Way sized galaxies. We find that cold mode accreted gas 
enters a galaxy halo with ~ 70% more specific angular momentum than dark matter averaged over 
cosmic time (though with a very large dispersion). In fact, we find that all matter has a higher spin 
parameter when measured at accretion than when averaged over the entire halo lifetime, and is well 
characterized by A ~ 0.1, at accretion. Combined with the fact that cold flow gas spends a relatively 
short time (1-2 dynamical times) in the halo before sinking to the center, this naturally explains why 
cold flow halo gas has a specific angular momentum much higher than that of the halo and often forms 
"cold flow disks" . We demonstrate that the higher angular momentum of cold flow gas is related to 
the fact that it tends to be accreted along filaments. 

Subject headings: galaxicslormation — halos — evolution mcthodsmumerical — hydrodynamic — 
simulation 



1. INTRODUCTION 

In the modern Lambda Cold Dark Matter (LCDM) 
paradigm of galaxy formation, dark matter halos form 
via the continuous accretion of diffuse material, as 
well as me rgers of sma l ler d a rk matter halos over 
time (e.g., iStewart et all 120081: iFakhouri fc"Mal 120091: 
IDiemer et al.l 120121 and references within). While 
the acquisition of angular momentum in dissipation- 
less iV-body simul a tions has been well s tudied (e.g. . 
Bullock et al.l I200TI: IVitvitska et al l [20021: IMaller et al 



2002t iD'Onghia fc Navarro! 120071 lAvila-Reese et al 



200a iBett et all 120071 120101 ), our understanding of 
the extent to which gas behaves in the same man- 
ner has o nly recently started to be explored compu- 
tationally (IChen et all 12001 ISh arma fc Steinmetz 2005 



Brook et al.ll20ToF IKimm et al. H2011l:IPichon et allboiTl 
Sharma fc Nathl I2012D In this paper, we extend such 



work by looking at angular momentum acquisition in 
galaxy halos in a suite of cosmological hydrodynamic 
simulations that allow us to directly measure the angu- 
lar momentum not only of dark matter, but also of the 
gaseous halos that are built from ongoing gas accretion 
from the cosmic web. 
In this LCDM framework, the canonical picture of 
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galaxy formation is that gas falling into a dark matter 
halo shock-heats to the virial temperature of the halo 
Silkl [19771 IWMFe fc Reesl [19781: IBarnes fc Efstathioul 



[19871 IWhite fc Frenkl I199l[ IMaller fc Bullockl I2004D 

spending a considerable amount of time in the galactic 
halo and redistributing angular momentum so that gas 
in the halo has the same specific angular momentum as 
the underlying dark matter. This shock-heated gas can 
then cool to the center of the halo where it presumably 
forms an angular momentum supported galactic disk. 
Under this formalism, extended gaseous halos around 
galaxies are expected to have similar specific angular 
momentum as the dark matter (by construction) and 
the size of galactic disks should be dir ectly proportional 
to the spirQ of the da rk matter halo ((Fall fc Efstathiou 
119801: IMo et al.l fl998l) . Indeed, this picture of gas ac- 
cretion and angular momentum acquisition is a vital 
component of most semi-analytic models of galaxy for- 
mation in order to model galaxies with realistic disk 
sizes based on underlying dark matter halo properties 
(recently, e.g.. iCattaneo et al.l|2006HCroton et alj|2006t 
iSomerville et all 120081: iDuttonl 120 121). However abun- 
dance matching techniques (e.g.. iBehroozi "eTall f2TTToh 
have demonstrated that at most ~ 20% of available 
baryons can be found in galaxies in the present day, 
thus the angular momentum of a galaxy may strongly 
depend on the origin of the baryons that form the disk 
([Mailer fc D ckcl 2002). In fact, recent simulations find 
that the gas in g alactic disks has higher spin than that of 
the dark matter (IChen et all 2003; S harma fc Steinmetzl 



2005tlBrooks et al.ll20l£ in agreement with the model of 



Mailer fe DekeiP f2002) 



One possible reason for this discrepancy is the treat- 
ment of gas accretion. Recent advances in hydrody- 
namic simulations as well as analytic arguments have 
stressed the importance of "cold-mode" accretion of gas 

9 For the sake of brevity, we often use the term "spin" throughout 
this paper as a synonym for specific angular momentum, except 
when specifically referring to the "spin parameter," A (see Eq. [TJ. 
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that does not shock-hcat to the virial temperature be- 
fore building the central galactic disk. This mode of 
accretion occurs when the cooling time of infalling gas 
is shorter than the compression time, making it diffi- 
cult to establish a stable shock for h alos masses below 
a critical valujB of ~ 1 12 M^ (e.g. iKeres et alll200l 
Dekel fe Birnboimll2006t iBrooks et al.ll2009t iDekel et al 



2009t IFaucher-Giguere fc Keresl 120101: IStewart et al 



201 lat Ivan de Voort et al.l I201U IFaucher-Giguere et al 
201 If ). This cold mode of accretion tends to have a 
very high angular momentum content JKeres et a,l.H2009t 
Keres fc Hernauist|l2009t lAgertz et all 200ft IBrooket al l 



2010tlPichon et al.ll201lt iTillson et al.ll2012D . with cold- 
mode gas in galaxy halos typically having ~ 2-5 
times more specific angular mome ntum than the corre- 
sponding dark ma tter in the halo ([Stewart c t alj|20~llbt 
IKimm et al]|20Tl . 

The goal of this paper is to investigate angular mo- 
mentum acquisition at moderate and low redshifts in a 
suite of fully cosmological high-resolution hydrodynamic 
simulations, in order to understand the origin of angular 
momentum in gaseous halos of galaxies — particularly for 
cold-mode gas — and to determine how angular momen- 
tum acquisition in these gaseous halos differs from the 
more thoroughly investigated spin of dark matter halos. 
While the build-up of high angular momentum gaseous 
halos invariably influences the evolution of the galaxies 
embedded at the center of each halo, we do not focus 
on the simulated disk galaxies' angular momentum in 
this work, as there is already a rich literature of low red- 



ilar to those presented here (e.g., 


Governato et al.ll200ft 


Brooks et al.l 12009 


; iBrook et al. 


20101: Stinson et al.l 


20101 IBrooks et al.l 


201lHGuedes et al.ll2011l). 



The outline of this paper is as follows. In section Sj2] we 
discuss the details of the cosmological hydrodynamic sim- 
ulations used and the effects that each simulated galaxy's 
unique merger history has on the more fundamental un- 
derlying trends of angular momentum acquisition. We 
analyze the angular momentum acquisition into galaxy 
halos over cosmic time in section $3j and present implica- 
tions for observations of extended gaseous disks of inflow- 
ing material in SjD We discuss our results in SJSJ including 
implications for implementing cool gas halos with realis- 
tic angular momentum content in semi-analytic models 
of galaxy formation. We conclude and summarize our 
findings in S|6] 

2. THE SIMULATIONS 

Each of the four simulations analyzed here uses a sepa- 
rate set of cosmological initial conditions to probe galax- 
ies in different environments, though all simulations use 
the same fundamental cosmological parameters, adopted 
from the bes t-fit parameters of the WMAP three-year 
data release (jSpergel et al.ll2007t) : M = 0.238, D, A = 
0.762, H = 73km s _1 Mpc -1 , n s = 0.951, and cr 8 = 
0.74. We focus our analysis on the single most massive 
galaxy in each simulation, tracking the evolution of each 
galaxy halo until z = 0. Each of our four galaxies reside 
in a roughly Milky Way size dark matter halo, and we re- 

10 "Cold-mode" accretion is also important in halos above this 
critical mass at high redshift, but we will focus on z < 2 for the 
majority of this paper. 



fcr to these halos as "Halo 1-4" throughout this papciO- 
Halos 1-2 both reach the critical halo mass for shock- 
heating infalling gas (M v ; r ~ 1O 12 M ) at some point dur- 
ing each simulation, with M v ii ~ 1.5 x 10 12 M Q at z = 
for both halos. In contrast, Halos 3-4 remain below this 
critical mass at all times, growing to M v j r ~8x 10 1:l A/q 
at z = 0. We define the virial radius of our halos as the 
radius within which the average density is A v j r times the 
mean density of the universe, iteratively removing un- 
bound particl es from within halos by run ning the Amiga 
Halo Finder (jKnollmann fc Knebd 120091) on all output 
snapshots of our simulatio ns^, where A v j r i s a fun ction 
of redshift and defined bv lBrvan fc Normanl (|1998l ). We 
note that this is a relatively standard method of defin- 
ing the virial rad ius of halos in dissi pationless A^-body 
simulations (e.g.. IStewart et aLl 120081 ). however there is 
also a growing concern that conventional definitions of 
halo virial radius and virial mass are somewhat arbitrary 
and that, ultimately, they may not be physically mean- 
i ngfu l (e.g. JAnderhalden fc Diemandll201ll : lDiemer et al.l 



l2012f ). As we are primarily concerned with comparing 
the angular momentum acquisition of galaxy halos via 
different modes (cold-mode gas, hot-mode gas, and dark 
matter accretion), the precise definition of the virial ra- 
dius used should not impact our results. 

We generate "zoom-in" multi-mass particle grid ini- 
tial conditions within a periodic box of 40(50) co-moving 
Mpc on each side for Halos 1-2 (3-4), in order to account 
for large scale tidal torques. The high-resolution parti- 
cle masses in the initial conditions are (md ar k, "?,gas) = 



(17,3.7) x 1O 5 M for Halos 1-2, and (m da 



k,m„ 







(12,2.1) x 1O 5 M for Halos 3-4. The simu lations use 
the code GASOLINE (|Wadslev et. all l200l . which is 
a smoothed particle hydrodynamics (SPH) extension of 
the pu r e A^-b ody gravity code PKDGRAV developed by 
iStadell (|2001l ). The gravitational force softening is 332 
(347) pc for Halos 1-2 (3-4), which evolves co- movingly 
until z = 9 and remains fixed from z = 9 to the present. 
The SPH smoothing length adapts to always enclose the 
32 nearest gas particles and has a minimum of 0.05 times 
the force softening length. 

The code assumes a unifor m UV backgroun d from 
QSOs, implemented following lHaardt fc Madaul ([1996D 
and F. Haardt (2002, private communication). It in- 
cludes star formati on, and Compton and radiative cool- 
ing as described in iKatz et al.1 (|1996| ). The "blastwave" 
supernova feedback model implemented creates turbu- 
lent motions in nearby gas particles, keeping them from 
cooling and forming stars; this model is implemented in 
the code by for cibly turning off cool ing for a short time, 
as described in iStinson et. all (|2006l ). The only free pa- 
rameters in this star formation and feedback model are 
the minimum density threshold (p m in = 0.1 cm -3 ), star 
formation efficiency factor (c* = 0.05), and fraction of 
supernova energy that couples with the ISM (esN = 0.1 
for Halos 1-2 and esN — 0.4 for Halos 3-4) which have 

11 Halo 1 uses the same initial conditions as the very high resolu- 
tion iV-body simulation Via Lactca II, whose merger history can be 
publicly downloaded at http : //www.uco lick . org/~diemand/v l/. 
Halos 1-2 are the same halos analyzed in Stewart ct al. (2011a b), 
while Halos 3-4 have also been analyzed in past works (referre d to 
as "h277" and "h285" respectively in, e.g., Brooks ct al. 2011) 

12 This unbinding decreases the virial radius and virial mass by 
only ~ 5%, compared to including unbound particles. 




Figure 1. Evolution of mass and angular momentum in the halos of 4 simulated galaxies over cosmic time for cold- mode gas (dotted 
blue line) hot-mode gas (dot-dashed red) stars (dashed green) and dark matter (solid black). Top: the total mass within the virial radius 
(after removing unbound particles). Bottom: specific angular momentum, in terms of the spin parameter, A (sec Eq. [TJ. Note how drastic 
changes in angular momentum often correspond to significant merger events (long vertical dashed lines in all panels) or even occasionally 
by gas-rich minor mergers (short vertical dashed line at z ~ 2 in Halo 2). 



all been motivated bv lGovernato et aTl ([20071 ) to produce 
galaxies with realistic star formation rates, disk thick- 
nesses, gas turbulence, and Schmidt law over a range in 
dynamic masses. Similar simulations utilizing this feed- 
back scheme have show n great success in prod ucing real- 
istic disk-type galaxies (IGovernato et al.l2009 ) , ma tching 
the mass-metallicity relation ([Brooks et al.l I2007T ) . and 
matching the abundance of D amped Lyman a systems 
at z = 3 (iPontzen et al.l [20081) . We refer the reader to 
iGovernato et al.l l 20091) for a more detailed description of 
the simulation code, though it is also worth noting that 
most of the analysis in this paper focuses on gas that has 
yet to form stars or be directly affected by supernova 
feedback. 

As we are focusing on angular momentum acquisition 
in this paper, we also note that simulations of this kind 
have produced galacti c disks of realistic s izes when com- 
pared to observations ([Brooks et al.ll20Tll ) . This is a par- 
ticularly important point, as it demonstrates that there 
is no significant artificial angular momentum loss that 
might lead to small disk sizes for our simulated galaxies. 

With the resolution implemented here, this feedback 
model results in minimal galactic winds (~ 100 km/s) 
that mostly affect hot gas, and are more pr ominent in 
halos with A/ vir < 10 n M o (jShen et all 120101 ) . Stronger 
outflows result when using high resolution simulations 
with higher star formation density thresholds and mod- 
ified star f ormation and superno v a efficiency param- 
eters (e.g.. IGovernato et al.l 120101: iGuedes et al.l [20111 : 
IGovernato et al.ll2012t iPontzen fc Governatol 120121 ). but 
the parameters adopted in this paper do not drive out- 
flows of cool gas. As a result, we may more easily focus 
our analysis on the behavior of cold-mode gas accretion, 
which is a fundamental first step to fully understanding 



the complex environment of gaseous halos around galax- 
ies, which are thought to include both accretion and out- 
flowing winds. 

In the analysis that follows, we define "cold-mode" gas 
particles based on a temperature cutoff of 250, 000 K, 
follow ing t he character i zation put forth by IKeres et al.l 
(|2005l ) and IKeres et al.l ([20091 ) . Gas that has never been 
above this temperature (until it reaches 0.1i? v ir) is de- 
fined to be "cold-mode." Any gas that does not meet 
this criterion is termed "hot-mode." We stop tracking 
temperatures at R = 0.1i? v j r so that these cold/hot mode 
definitions are not contaminated by heating due to super- 
nova feedback. We do not make any special distinction 
for cold-mode gas that is associated with mergers. 

We note that smooth particle hydrodynamic (SPH) 
codes are known to have some nume rical issues fo l lowin g 
cool clouds in a hot medium (e.g., lAgertz et al.l I2009T ) . 
however most of the general results about cold mode ac- 
cretion have also been confirmed in grid-based adapted 
mesh refinement (AMR) simulations, though many small 
scale details arc different. We believe the results in this 
paper in regards to angular momentum are robust to nu- 
merical implementation as we are mostly concerned with 
the state of gas in the the outer halo where resolving the 
multi-phase medium is less of a concern. In addition, 
iPichon et al.1 ([2011D have per formed a simila r study us- 
ing the AMR code Ramses ([Tevssierl I2002D , and their 
results are consistent with ours. 

2.1. Mergers and Stochasticity 

The top panels of Figure [T] show the mass growth as a 
function of time for the dark matter as well as the various 
baryonic components within the halo (all material within 
the virial radius, including the galaxy itself). We note 
that there is significantly more cold-mode gas in galaxy 




Figure 2. Specific angular momentum of accreted material. The dashed green curve represents a semi-analytic model where halo virial 
quantities are known, and a constant spin parameter is assumed, such that j(t) = A\/2V^ir(t)-Rvir(i)- Top panels show only fresh accretion — 
material that entered the virial radius within ~ 300 Myr. We see that cold- mode gas on average has ~ 50% higher spin than all gas at 
th e time of accretion. We compare our fixed spin parameter model (A = 0.1, green clashed curve) to the more complex tidal torque model 
of IMaller et al.l Q2002D . given by the gold dashed curve. The bottom panels look at all material currently in the halo, compared to the 
canonical model where A = 0.04. 



halos than there is gas in the galaxies themselves, for all 
four simulations and at all epochs. Thus the cold-mode 
gas can largely be thought of as warm halo gas, like tha t 
seen in QSO absorption systems (e.g.. lWerk et al.ll2012f ). 

The bottom panels of Figure Q] show the specific an- 
gular momentum of gas and dark matter in each galaxy 
halo. Unlike the top panels, we now remove the galaxy 
itself from consideration, and only look at material with 
0.1i?vir < R < Rvh- We quantify the specific angular 
momentum in this figure in terms of a dim e nsionl ess spin 
parameter, A, adopted from lBullock et ail (|2001[ ). as: 



where Aa; is the spin parameter of a particular compo- 
nent x (dark matter, hot-mode gas, or cold- mode gas), 
j x is the specific angular momentum (J/M) of that com- 
ponent within a sphere of radius R, and V is the circular 
velocity measured at R. As in past works, we again find 
that the spin of cold-mode gas in the ha lo is typically 
much higher than that of the dark matter (jStewart et al.l 
l2011btlKm7m et al.ll27)lTh . 

Throughout this paper, we calculate the specific angu- 
lar momentum of each component (dark matter, cold- 
mode gas, hot-mode gas) in relation to the spin axis 
of that same component. In this way, we may more 
robustly compare the angular momentum of gas versus 
dark matter without artificially reducing the spin of one 
component by projecting it along the axis of a different 
component. This approach also makes the most sense 
observationally, as it is not feasible to project any ob- 
served angular momentum measures of gas in the ha- 
los of galaxies with the spin axis of the underlying dark 



matter, which cannot be observed. Thus, when the cold- 
mode and hot-mode gas in our simulations show similar 
spin parameters in the bottom left panel of Figure Q] (for 
Halo 1 at z < 1), this does not necessarily mandate that 
each gaseous component is spinning along the same axis 
at this time, but that the total amplitude of each compo- 
nent's J/M is similar, regardless of the relative alignment 
with respect to each other (we explore the non-trivial re- 
lation between the spin axis of these various compo nents 
in an upcoming paper, llvloustakas fc Stewart! [20ll . 

The tall vertical dotted lines in these figures show when 
each galaxy experiences a significant merger eveniP^I and 
often accompany sharp peaks in the mass growth of the 
galaxies. Often, these mergers significantly re-orient the 
angular momentum axis of the galaxy, the halo, or both, 
indicated by large peaks or dips in the bottom panels of 
Figure [TJ The short vertical dotted line in Halo 2 shows 
a minor gas-rich merger that seems to have an unusually 
strong effect on the angular momentum of the gaseous 
halo. It is clear in the bottom panels that the spin of 
the halos is quite stochastic — especially for the cold gas 
in the halo. In a more detailed sense, these panels show 
that the impact a merger has on the angular momentum 
of the galaxy halo depends on the precise characteristics 
and orbital parameters of the merger. There are multi- 
ple instances where mergers result in a drastic increase 
of the specific angular momentum of material in the halo 
(for example, Halo 1 and Halo 4). However, there are 
also times where mergers result in a rapid decrease in 

13 A merger of stelar mass ratio greater than 0.3. These lines 
correspond to the time of final coalescence, with the understanding 
that mergers are not instantaneous events. It is likely that each 
merger impacts the halo for several billion years prior to the times 
indicated in Figure [Tl 
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the angular momentum of the halo, due to misalignment 
between the pre-existing spin axis of the halo and the 
orbital parameters of the incoming merger (Halo 2 at 
z ~ 0.8 and Halo 3). While the detailed impact that 
these merger events have on the angular momentum of 
the halo would be an interesting topic of future study, 
the goal of this paper is to understand the underlying 
steady acquisition of angular momentum into galaxy ha- 
los over cosmic time, apart from the stochastic peaks and 
dips that result from merger events. As a result, we will 
not be focusing on the detailed minutiae of the angular 
momentum growth in each of our galaxy halos, which 
will depend strongly on the individual accretion history 
of each halo. 

3. ANGULAR MOMENTUM ACQUISITION 

We first look at how the angular momentum of fresh 
gas accretion to the virial radius compares to the canon- 
ical scenario for dark matter. In the top panels of Fig- 
ure [2] we directly calculate the specific angular momen- 
tum of recent accretion to the virial radius, for both gas 
and dark matter as a function of cosmic time, defining 
"recent" accretion by a timescalc of /lesssimSOO Myr, 
which is shorter than the halo dynamical timescale at all 
epochs of interest (z < 3; see ^3. II for more on gas accre- 
tion timescales in the simulations), but long enough to 
be practical in terms of storing simulation outputs with 
a finite timestep resolution. Particles are considered "ac- 
creted" onto the halo once they enter the virial radius of 
the halo. On occasion, a particle may enter, exit, and 
re-enter the halo at a later time; in this scenario, only 
the latest entry to the halo is considered to be true ac- 
cretion to the halo, to avoid double-counting the angu- 
lar momentum contribution of any given particlq^j. For 
comparison, the lower panels of Figure [5] show the to- 
tal specific angular momentum of all material within the 
virial radius, including the galaxy itself. 

It is well documented that dark matter halos in N- 
body simulations show r elatively stable spin parame- 
ters of A ~ 0.04 (recently. IBett et al.ll2007t IMaccio et all 
[20071: IBerta et al.|[200l IBett et alJl2uToir Since the spin 
parameter is defined to vary as X(t) oc j(t)V^ T i2~ ir , and 
since we know that for any given halo the circular velocity 
(Vvir) and the virial radius (i? V ir) both grow as a func- 
tion of cosmic time, the only possible way to maintain 
a constant spin parameter is if the specific angular mo- 
mentum, j(t), increases over cosmic time. As a result, it 
is a robust prediction of LCDM that the specific angular 
momentum of freshly accreted material must, on aver- 
age, be higher than th a t of the halo as a whole (also see 
Porciani et all |2002albt IMaller et aJ]|2002t iPichon et all 
2011ft . This is an important context to keep in mind for 
much of the following analysis. 

In each of the top panels, we compare our findings to a 
simple model where a constant spin parameter of A = 0.1 
is assumed (i.e., j ~ 0.14i? vir V V ir) given by the dashed 
(green) curve, noting that this is ~ 2.5 times higher 
than the global spin parameter of the dark matter. We 
also com pare to the more co mplex model (gold dashed 
curve) of IMaller et al. (2002j), in which the angular mo- 
mentum of fresh accretion, in radial shells around the 

14 At any given time, < 20% of all particles that ever entered 
R v i r have exited the virial radius of the system. 
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Figure 3. Specific angular momentum distribution. P(j) of dark 
matter, cold-mode gas, and hot-mode gas at R = R V1L - around Halo 
2 at z ~ 1.4. To aid in comparison between components, the x-axis 
has been normalized by the average specific angular momentum of 
the dark matter at the virial radius 

halo, has been acquired at t urnaround according to lin- 
ear tidal torque th eory (e.g., Peebles 19691 iDoroshkevichl 
Il970t IWhitelll984[ ). Remarkably, the simple assumption 
of A 0.1 for infalling material produces very similar 
curves to the tidal torque model, reproducing the spe- 
cific angular momentum of fresh accretion as well as (or 
better than) the physically motivated tidal torque model. 

In the bottom panels, the expectation of A = 0.04 is 
a much closer approximation to the curves. The differ- 
ence between the top and bottom panels of this figure 
(a factor of ~ 2.5 in j) is paramount for understanding 
the angular momentum content of galaxy halos versus 
that of galaxies. We see from this figure that the an- 
gular momentum of all baryons is very similar to that 
of all dark matter, both in fresh accretion (top) and in 
all material within the virial radius (bottom); thus, the 
canonical assumption of A = 0.04 still seems to be a rea- 
sonable estimate for the angular momentum of galaxies, 
or of all baryons in the entire halo. However, this as- 
sumption grossly underestimates the spin of recent dark 
matter accretion into the halo, as well as that of the 
gaseous component in the halo. Note that for each of 
our simulated galaxies, both in fresh accretion as well 
as cumulative accretion, cold-mode gas is almost always 
higher than that of recent dark matter or even that of 
hot-mode gas accretiorP^I. The cold-mode gas has on av- 
erage a ~ 70% higher spin than that of the dark matter at 
the time of accretion, though with very large variations. 
This systematic offset between the angular momentum 
of cold-mode versus hot-mode gas accretion suggests a 
fundamental kinematic difference between gas that is ac- 

15 A notable exception to this result is for Halo 3 at z ~ 1.5, 
when j for cold-mode gas accretion is less than that of the recent 
dark matter accretion, presumably as a result of the mergers in 
this systems shown in Figurefl] While some authors distinguish be- 
tween cold-mode gas accretion and cold gas associated with merger, 
we remind the reader that we make no such distinction here. 
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Figure 4. Example velocity flow chart for recent gas accretion. In each panel, the length of each mark is proportional to the projected 
velocity of the corresponding particle, while the colors correspond to different accretion times to the virial radius of the halo. The width of 
each panel is ~ 200 co- moving kpc (the virial radius for this halo at this epoch), with left and right panels showing orthogonal projections 
of the same galaxy. The galaxy is viewed edge-on in the left panels and face-on in the right panels, while the top panels show cold-mode 
accretion and the bottom panels show hot-mode accretion. 



creted in each of these modes. 

To gain further insight into how these modes of gas 
accretion differ from each other, and from dark matter 
accretion, Figure [3] shows a representative mass- weighted 
histogram of the specific angular momentum for particles 
of dark matter (solid black), cold- mode gas (dotted blue), 
and hot-mode gas (dot-dash red), focusing on a single 
simulation at a single point in time (Halo 2 at z ~ 1.4). 
For the sake of straightforward comparison, we normal- 
ize the x-axis by the average specific angular momentum 



of all dark matter at i? v ir, <3dm (Rvit)>, and the y-axis 
by the cosmic baryon fraction. We show the distribution 
P(j) for particles between 0.95i? vir < R < 1.05i2 v i r . The 
vertical line shows the mean value of P(j). We see that 
the specific angular momentum distribution of dark mat- 
ter is much broader then that of the gas components. The 
hot-mode gas is much narrower and has a mean near zero, 
which we might expect if the hot gas is shock heating as 
it enters the virial radius. What is more surprising is that 
the cold-mode gas has a distribution very different than 
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the dark matter with a much higher mean, even though 
we do not expect the cold-mode gas to have exchanged 
angular momentum in any coherent fashion. This dual- 
peak distribution of cold-mode gas at this epoch is not 
due to a gas-rich merger or some other unusual accre- 
tion event, but seems to be a typical result of coherent 
cold-mode gas accretion along cosmic filaments. We have 
checked that these differences remain at 2i? vu - suggest- 
ing that they have to do with the large scale distribution 
of the various components and not particle interactions. 
Taken alongside Figure [2j this suggests a fundamental 
kinematic difference between the origin of cold-mode gas 
versus hot-mode gas accretion onto galaxy halos. 

We investigate this distinction further by visualizing 
all recently accreted gas particles within the virial radius 
of this same halo at this same epoch in Figure |3J Each 
panel shows the position of all cold-mode (top) and hot- 
mode (bottom) gas particles that were accreted to the 
halo within the past 2 Gyr, with the length of each mark 
proportional to the projected velocity of each particle, 
so that the figure gives a qualitative visual impression 
of the velocity flow of the particles shown. The color of 
each particle corresponds to the accretion time, allowing 
a sense of time progression from the most recent accretion 
(blue; T < 300 Myr) to the older accretion that has had 
time to sink to the center of the halo and mix with the 
pre-existing gaseous halo (red; 1 < T < 2 Gyr). The 
galactic disk is edge-on at the center of the image in the 
left panels, and face-on in the right panels, as indicated 
by the red particles in the top plots where cold-mode 
gas has had sufficient time to reach the galaxy. In the 
bottom panels, the galaxy cannot be seen, as hot-mode 
gas accretion typically takes longer than 2 Gyr to reach 
the galaxy at this epoch (we discuss accretion timescales 
further in g3T|l . 

We note several important qualitative features in Fig- 
ure [5] The most striking is the stark difference between 
cold-mode and hot-mode accretion. As has been noted 
in past simulation studies, cold-mode accretion typically 
flows into the halo along the direction of the large-scale 
filaments that fuel the galaxy, most clearly visible in the 
upper-left panel, while ho t-mode accretion i s more likely 
to be more isotropic (e.g. jKeres et al.ll2009f) . The veloc- 
ity vectors of the particles also show that while both cold- 
mode and hot-mode gas accretion has significant angular 
momentum on a particle-by-particle basis, the cold-mode 
gas in the upper panels shows clear indication of coherent 
rotation. While no such bulk motions are readily visible 
by simple visual inspection in the lower panels, Figure [5] 
has already demonstrated that the hot-mode gas is also 
spinning about the galaxy with significant specific an- 
gular momentum — if not quite as high as is present in 
cold-mode material. 

The upper panels also show clear segregation of cold- 
mode gas particles as a function of accretion time. As 
gas first accretes onto the halo, it clearly begins in the 
outskirts near the virial radius (blue) and gradually falls 
inward, spinning about the galaxy and losing orbital en- 
ergy as it sinks towards the center of the halo (green- 
gold). After 1-2 Gyr, the cold-mode gas reaches the 
galactic region and begins to fuel the outer edges of the 
galaxy itself (red). This clear picture of gas accretion 
from the cosmic web directly onto the galactic disk is 
not indicated by the bottom panels for hot-mode accre- 
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Figure 5. Spin parameter of gas and dark matter, revisited. The 
colored lines show the spin parameter of recent accretion (T < 
0.3 Gyr, R < 0.5i? v i r ). When only recent dark matter accretion 
along cosmic filaments is included, it has a much higher angular 
momentum than the fiducial value of A ~ 0.04, and is more similar 
to the spin of recent cold-mode gas accretion, which also flows into 
the halo along filaments. 

tion. Instead, we see significant radial mixing, with gas 
from a wide range of accretion times at various radii. 

If the different spatial and kinematic origins of filamen- 
tary (cold-mode) accretion versus non-filamentary (hot- 
mode) accretion is truly what leads to the distinct angu- 
lar momentum distributions on infall, then we should see 
a similar distinction between the spin of filamentary ver- 
sus non-filamentary dark matter accretion as well. Fig- 
ure [S] shows the spin parameter as a function of time 
for recently accreted (T < 300 Myr) dark matter, cold- 
mode gas, and hot- mode gas, using an identical color/line 
scheme as in Figure [1] However, we have now broken up 
the recent dark matter accretion into two mutually ex- 
clusive components: filamentary dark matter accretiorPI 
(black solid curve) versus non-filamentary dark matter 
accretion (black dashed curve). Aside from some no- 
ticeable stochasticity, we see remarkable correlation be- 
tween the specific angular momentum of fresh filamen- 
tary accretion (filamentary dark matter versus cold-mode 
gas) and between that of fresh non-filamentary accre- 
tion (non-filamentary dark matter versus hot- mode gas) . 
These correlations are highly indicative that the higher 
specific angular momentum content of fresh cold-mode 
accretion is a fundamental consequen ce of filamentary ac - 
cretion to dark matter halos (see also lPichon et al.ll2011l) . 

Because filamentary accretion is often thought to nec- 
essarily indicate highly radial infall trajectories, we note 

16 At each epoch, we define filamentary accretion spatially by 
locating particles in a bi-conical region, extending from the halo 
center towards the direction of the filament (determined by visual 
inspection of matter density contours) in either direction, with a 
3D opening angle of 20 degrees. 
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Figure 6. The average amount of time it took for gas currently in the galactic region (_R < 0.1-R v i r ) to sink from -R v i r to 0.1R, v i r . Left: 
Note that the sinking time for cold-mode gas accretion is relatively short, with timescales typically between 1 — 2r, where r is the halo 
dynamical times and is given by the two solid green curves in both panels. Right: The sinking time for hot-mode gas accretion is typically 
~ 1 — 2 times longer than it is for cold-mode gas. The median lookback time to accretion for dark matter particles in the halo is considerably 
longer (~ 7 — 10 Gyr at z = 0). 



that this kinematic difference between hot-mode and 
cold-mode gas is not simply due to their average orbital 
parameters upon infall. Quantifying the "typical" orbit 
of fresh accretion to the virial radius in terms of a median 
infall angle — the angle between a completely radial in- 
fall vector and the actual velocity vector of each particle 
upon first entering the virial radius — we find that cold- 
mode gas, hot-mode gas, and dark matter all have com- 
parable mean infall angles ranging from ~ 30-50 degrees 
from radial, with no clear trends for any component being 
consistently more or less radial than other components. 
We want to take a moment to emphasize this point, as 
it has important ramifications for absorption-line stud- 
ies of gaseous halos. When interpreting velocity offsets 
in studies of absorption systems along lines of sight to 
background galaxies, it is often assumed that infalling 
gas accretion onto the galaxy should always be indicated 
by high-velocity pure radial infall. We emphasize that 
this is not the case; for both dark matter and gas, the 
vast majority of accreted material enters the virial ra- 
dius with high angular momentum, spinning about the 
halo at an initial angle of ~ 40 degrees from a purely ra- 
dial trajectory, with only a portion of the infall velocity 
along the radial direction. As a result, we advise that 
absorption-line searches for cold-mode accretion should 
not focus exclusively on signs of purely radial infall, but 
should also look for signs of coherent rotation in cool 
gaseous halos (for obse rvable signatures of c oherent ro- 
tation in Halos 1-2, see lStewart et al.ll2011b[ ). 



3.1. Halo Gas Accretion Timescales 

Figure 2] indicated that the radius of infalling cold- 
mode material is tightly coupled with its accretion time, 
whereas the picture is more complicated for hot-mode gas 
(or dark matter accretion). In order to understand the 
processes involved in transporting high angular momen- 
tum fresh accretion from the virial radius to the galaxy, 
it is important that we recognize the different timescales 
involved in this transport. However, we first must ask the 
question: once infalling material falls to a given radius, 
does it stay within that radius? If material routinely en- 
ters the galactic region only to return along an elongated 
orbit to the outer edges of the halo again, it will be dif- 
ficult to define an "infall timescale" in any meaningful 
way. 

We compare the cumulative mass of all material that 
was ever accreted to within a given radius, R, to the mass 
currently enclosed within R, and find that the vast ma- 
jority of gas and dark matter that falls within the virial 
radius is ultimately bound to the halo, with < 20% of all 
mass that ever entered the halo being lost to the system, 
most likely during violent major mergers. In contrast, 
when we focus on material that has fallen within the in- 
ner halo (R < 0.5i? v ir) we notice an important difference. 
While ~ 85-95% of all gas that has ever entered the in- 
ner halo remains there, only ~ 60-70% of dark matter is 
equally bound to this inner radius. 

The difference between this behavior of dark matter 
and gas within the halo is not surprising, as the highly 
elliptical elongated orbits of collisionless dark matter par- 
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tides will often bring them back to the outer halo for 
significant periods of time during each orbit about the 
halo's center. In contrast, high angular momentum gas 
particles are more likely to mix with other gas in the in- 
ner halo, redistributing angular momentum so that they 
no longer have enough orbital energy to exit the inner 
halo. Consequently, while it is indeed straightforward to 
define a "sinking time" for inf ailing gas, it is not mean- 
ingful to define a similar quantity for dark matter. Thus, 
for each halo we define the sinking time by tracking every 
gas particle currently in the galactic region of the halo 
at each epoch (R < 0.1i? V ir) back in time in the simu- 
lations, calculating the amount of time it took to reach 
R = 0.1i? v i r after first entering the virial radius. 

Using this definition, Figurc[6]shows the mass- weighted 
mean sinking time for cold-mode gas particles in each 
of our simulations as a function of time. Because the 
sinking timescale invariably depends on halo properties, 
which in turn depend on the epoch in question, we find 
it useful to compare these timescales to the halo dynam- 
ical timf3r(z). We compare the sinking time to r and 
2r (the solid green curves), finding that typical sinking 
times are relatively short for cold-mode gas — only ~ l~2r 
(< 3 Gyr at z = 0). We repeat this analysis to deter- 
mine the average sinking time from R = 1.0-0.5i? v h, for 
all material currently in the inner halo (R < 0.5i? v ir)- 
Not surprisingly, the sinking time for material still in the 
inner halo is roughly half of the values shown in Figure O 
These findings strongly suggest that infalling cold-mode 
gas in galaxy halos is not a cumulative reservoir of an- 
cient past accretion, but instead probes fresh accretion 
to the halo, with very high specific angular momentum 
content similar to that presented in the top panels of 
Figure [2] 

The sinking time for hot-mode gas accretion is typ- 
ically longer than that of the cold-mode curves shown 
here by a factor of ~ 1 — 2. As mentioned previously, 
we cannot provide a true "sinking time" comparison for 
dark matter particles, but we instead calculate the me- 
dian look-back time to first accretion for dark matter par- 
ticles currently in the halo at z = and find considerably 
longer timescales: ^7—10 Gyr for dark matter particles 
in 0.1R vir < R < 1.0-Rvir and ~ 5 — 8 Gyr for dark matter 
particles in the outer halo (0.5i? V h- < R < 1.0i? v ir)- 

4. OBSERVATIONAL CONSEQUENCES 

Based on our findings, we know that fresh cold-mode 
gas accretion onto galaxy halos has much more specific 
angular momentum than is typical for the entire dark 
matter halo (Figure [H top), or for all baryons in galaxies 
(Figure [5J bottom) , and that fresh cold- mode gas accre- 
tion even has higher specific angular momentum than 
recent hot-mode gas accretion, due to its fundamentally 
different origin along cosmic filaments (Figures EH5]). It 
is possible that these inherent difference in the angular 
momentum content of cold-mode versus hot-mode accre- 
tion might play a role in explaining why spiral and el- 
liptical galaxies appear to retain different percentages of 

17 We define the halo dynamical time by the crossing time: r = 
R/V oc (A„(z) p u (z))~ - 1 / 2 , such that t is independent of halo mass 
and only depends on reds hift. We adopt simpli fied fitting functions 
for t(z) as presented in IStewart et all 1)20091) . At redshifts z = 
3, 2, 1, 0, t(z) ~ 0.3, 0.5, 1, 2 Gyr respectively. 



the specific angular mome ntum of all infalling material 
(|Romanowskv fc Fall|[2012l ). 

In addition, we find that infalling gas does not travel 
through the halo in purely radial orbits, but rather has 
typical infall angles of ~ 40 degrees (with respect to ra- 
dial), with the spin of cold- mode gas in the halo often in a 
coherent direction, set by the direction of the large-scale 
cosmic filaments (Figures [4]{5j) . We have also established 
that cold-mode gas in galaxy halos has relatively short 
sinking timescales, and thus probes fresh accretion (Fig- 
ure |6]). 

Taking these results together, it is not surprising that 
cool halo gas around galaxies often forms coherent struc- 
tures with very high angular momentum, formed from 
continually infalling material from the cosmic web that 
eventually fuels the galaxy. As a result, we find that 
cold-mode gas accretion often forms extended "cold flow 
disks" of dense gas that rotates coherently about the 
galaxy, in roughly the same direction as the galaxy. Pre- 
vious simulation studies have demonstrated that this gas 
accretion should be observable as extended, warped disks 
of cool gas dense enough to form neutra l hydrogen (e.g., 
iRoskar et al.l2010t IStewart et alJl2011b|). which may help 



expla in observations of XUV d isks dThilker et"aE 



20071) and extended HI disks (E arcia-Ruiz et al 



2005 



2002 



Oosterloo et~aLll2007t iWalter et al.l 120081 ) around galax- 



ies, which tend to co-rotate with the inner galactic diks 
(|Christlein fe Zaritsky|[200l . 

While cold flow disks are not universal phenomena at 
all times in our simulations, neither are they rare iso- 
lated incidents. Figure [7] gives a visual impression of the 
rotation of cool gas within R < 100 co-moving kpc over 
a broad redshift range and sampling each of our four 
simulations, instead of only focusing on a single halo at 
a single epoch as we did in Figure SJ For each pixel 
in the images, the blue-red shading corresponds to the 
mass-weighted average velocity of HI gas along the line 
of sight, where we only show sightlines above a mini- 
mum detectability threshold of A/hi > 10 16 cm" 1 . All 
four galaxies show clear signs of coherent rotation of gas 
in the halo, as well as co-rotation between the gaseous 
halo and the galactic disk — even out to large radii. The 
co-rotation of the galaxies and their cold flow disks are a 
natural consequence of the high angular momentum con- 
tent of cold-mode gas in the halo, coupled with the fact 
that cold-mode gas accretion quickly sinks to the center 
of the halo, feeds the galactic disk, and thus helps define 
the angular momentum direction of the galaxy over time. 

Based on our simulations, we predict that cold flow 
disks (possibly already observed as XUV disks or ex- 
tended HI disks) should be more common for galaxies be- 
low the critical shock- heating threshold M v ; r < 1O 12 M , 
though they will likely extend farther into the halo 
for ga laxies very close to this threshold (|Stewart et al.l 
120 Hal) . They should contain more specific angular mo- 
mentum than their associated galaxies or dark matter 
halos, and are likely to be more metal-poor than the as- 
sociated galaxies (since the material is fresh accretion 
from the cosmic web). In addition, since cold flow disks 
probe recent cold-mode gas accretion to the halo, the 
orientation of these disks (often misaligned with the cen- 
tral galaxy) may give better indication of the direction 
of large scale cosmic filaments in which halos are embed- 




Figure 7. "Cold flow disks" in each of our four simulated halos, at 2 ~ 1.6, 1.4, 0.7, 0.2 respectively going from the top left to the bottom 
right panel. All panels are on the same velocity and co-moving distance scale, with the width of each panel set to 200 co-moving kpc. Note 
the coherent rotation in the cool halo gas well beyond the inner regions where the galaxy resides. The galactic disk is viewed near edge-on 
in each panel (i> 65). 



ded than the orientations of galaxies themselves We will 
expand on this conjecture fu rther in an upcoming paper, 
iMoustakas fc Stewart! (|2013l ). 

5. DISCUSSION 

In the canonical picture of angular momentum acqui- 
sition of galaxies, it is assumed that the angular mo- 



mentum of baryons (galaxies) are set by the underly- 
ing spin of t he dark matter halo in which each galaxy 
resides (e.g., IMo et ail 119981) . This formalism has led 
to great successes in reproducing galaxy population us- 
ing semi-analytic models of galaxy formation in which 
galaxies and their properties arc added on top of under- 



lying A^-body simulati ons of dark matter structure (e.g., 
ISomerville etall[200l . 

The assumptions inherent in this canonical picture of 
angular momentum acquisition do not work for gaseous 
halos of galaxies. While galaxies represent a cumulative 
growth of material over long time scales (similar to the 
underlying dark matter), gaseous halos probe the most 
recent accretion of fresh material from the cosmic web — 
especially for cool gas resulting from cold-mode accre- 
tion, which has very short sinking times of only a few 
billion years. To further complicate matters, cold-mode 
and hot-mode accretion contribute different amounts of 
specific angular momentum, with cold-mode accretion 
having higher specific angular momentum upon infall due 
to its distinct filamentary nature. 

Based on the results presented here, we propose a new 
scenario for angular momentum acquisition, in which: 
1) cold- mode gas typically carries more specific angular 
momentum than hot-mode gas accretion — a fundamen- 
tal difference that appears to be linked to its filamentary 
accretion nature 2) the sinking time of cold-mode gas 
is quite short, comparable to the halo dynamical time 
3) fresh accretion from the cosmic web (and thus, gas 
present in the halos of galaxies) has a much higher spin 
than that of the galaxy or the dark matter halo, spin- 
ning about the halo as it falls in, often in the form of 
extended cold flow disks 4) it is this high angular mo- 
mentum gas that provides galaxies with both mass and 
angular momentum, helping to set the spin axis of the 
galaxy. 

We note that in detail the full picture of angular 
momentum acquisition will be more complex than the 
scenario outlined above, in part due to the presence 
of galaxy scale gaseous outflows, which have been em- 
phasized in recent observations of cool halo gas (e.g., 
iSteidel et al1ll996l : lMartir][200l [Rubin et al.ll2010h and 
were not included in the simulations presented in this 
paper. It is worth noting however, that a recent simula- 
tion bv lBrook et~aLl (|201O) — one that includes bi-conical 
galaxy-scale outflows as a result of their high resolution 
and modified feedback scheme — found that the high an- 
gular momentum inflowing cold-mode gas was least af- 
fected by the outflowing gas. This was primarily due 
to the fact that inflowing material was more likely to 
approach the galaxy along its major axis, while out- 
flowing material was more likely to contain low angu- 
lar momentum gas ejected along the minor axis of the 
galaxy. Thus, the two mechanisms — cold-mode inflow 
and ejected outflows — typically inhabit different spatial 
regimes within the halo. 

5.1. Serai-analytic Models of Galaxy Formation 

While recent semi-analytic models have begun im- 
plementing co ld-mode and hot-mode gas accretion into 
galaxies (e.g.. iCroton et al.l l2006t iCattaneo et all 120061: 
ISomerville et all 12008ft . most of these models still pre- 
sume that the angular momentum of galaxies are pri- 
marily set by the spin parameter of the dark matter 
within the halo. Though this may provide a reason- 
able approach for estimating stellar disk sizes (because 
these are the baryons that were accreted early enough 
to have already formed stars), it will be difficult (if not 
impossible) to explain the presence of massive extended 
HI disks around galaxies with typical halo spin pa- 
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rameters (iGarcia- Ruiz et al.|[2002l: lOosterloo et al.ll2007t 
iWalter et al.l 12008ft or to accurately model the velocity 
structure of cool gas halos in the circum-galactic regions 
of galaxies without accounting for the fact that recent 
accretion is spinning more quickly than the bulk of the 
dark matter in the halo, 
i 

We propose that semi-analytic models should consider 
estimating the angular momentum of cool gas halos rel- 
atively simply by adopting the higher spin parameter 
(some ~ 3 times higher than the canonical spin of the 
halo) for all recently accreted material, motivated by the 
top panels of Figure 

Since this cold-mode halo gas is continually accret- 
ing onto the galaxy itself, one could also model the an- 
gular momentum of fresh gas accretion onto the outer 
edges of the galaxy (that is, cold-mode gas that is cur- 
rently falling on to the galaxy for the first time) by an 
angular momentum j(t) ~ X\/2V vlT (t — T)R v i r (t — t), 
where A ~ 0.1, and r is the sinking time (Figure 
comparable to the halo dynamical time. However, we 
note that it is non-trivial to integrate this continuous 
fresh accretion (to the galactic region) to obtain the to- 
tal angular momentum of the resulting galaxy. Over 
time, this high-spin material experiences significant vec- 
tor cancelation, due to the change in relative orienta- 
tions between the galaxy, the dark matter, the cold- 
mode halo gas, and the large-scale cosmic filaments 
fro m which new mate rial is constantly being accreted 
(seeE imm et all 120111 ). These orientations do not stay 
fixed over cosmic time, but vary somewhat stochastically 
based on large scale tidal fields and non-linear torques 
during major mergers. Instead, the direction between 
the galaxy orientation, the entire halo, and that of the 
inflowing cosmic filaments are more likely to be offset by 
~ 20 — 40 degrees, instead o f being perfectly aligned (e.g. , 
van den Bosch et all l2002t iSharma fe Steinmeti 120051 : 
Kimm et al. l l2011HMoustakas fc Stewartil2013ft . 

6. CONCLUSION 

In this paper, we used four high resolution cosmologi- 
cal hydrodynamic simulations to investigate the acquisi- 
tion of angular momentum in roughly Milky Way sized 
galaxy halos over time. Motivated by recent findings 
that the specific angular momentum of gaseous halos 
are c onsiderably higher than that of the da rk matter 
halo (|Stewart et al.l 1201 lbb iKimm et all 120 lift we focus 
on how reservoirs of cold-mode gas in galaxy halos grow 
and evolve over time, and how the angular momentum of 
this material differs from the canonical scenario for an- 
gular momentum growth in galaxies. We note that the 
angular momentum history of galaxy halos is sometimes 
quite stochastic — especially during galaxy mergers that 
may reorient the angular momentum axis of the system — 
so we focus on the average behavior of each system over 
time. We summarize our findings below. 

1. Cold- mode gas enters the halo preferentially along 
the direction of cosmic filaments, while hot-mode 
accretion is more isotropic in origin. As a result, 
cold-mode gas inherently has more specific angular 
momentum at first infall than hot-mode gas accre- 
tion (~ 70% more) . When filamentary dark matter 
accretion is compared to non-filamentary accretion, 
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we find similar differences between the initial spe- 
cific angular momentum of dark matter upon infall. 

2. The spin of all gas within the halo (including the 
galaxy) is broadly consistent with well-studied dark 
matter simulations, which find a relatively stable 
value of A ~ 0.04 over a broad range in epoch 
and mass. However, when we focus only on fresh 
accretion — gas and dark matter that has only been 
in the halo for < 300 Myr — we find much higher 
spins, consistent with A ~ 0.1. 

3. While gas accretion onto galaxy halos is often pic- 
tured as purely radial infall, we find that this is 
universally not the case in our simulations. Due to 
the high angular momentum of infalling material, 
it always spins about the halo with a typical infall 
angles (with respect to radial) of ~ 30 — 50 degrees 
at the virial radius. 

4. While the elongated orbits of collisionless dark 
matter particles may allow dark matter to popu- 
late the halo (even at the outskirts, near the virial 
radius) for long periods of time, inflowing gas is rel- 
atively short-lived in the halo. The "sinking time" 
for cold-mode gas to travel from the virial radius 
to accretion to the galactic region (R < 0.1i? V ir) is 
~ 1-2 halo dynamical times (< 3 Gyr at 2 = 0). 

5. This high angular momentum, quickly sinking, fila- 
mentary cold-mode accretion often takes the form 
of extended cold flow disks of newly accreted in- 
flowing material, coherently rotating roughly in 
the same direction of the galaxy. These extended 
gaseous disks may be related to observations of 
XUV and extended HI disks around galaxies. 

The nature of gas cooling in simulations invariably de- 
pends on the exact simulation code used, with different 
codes producing slightly different results from one an- 
other depending on the implementation and resolution, 
as well as different schemes for modeling the effects of 
supernova feedback, and in what way it injects energy 
back into the ISM. As a result, the definition of what gas 
is considered cold-mode versus hot-mode could be ar- 
gued to be a somewhat arbitrary distinction, especially 
when defined by a single maximum temperature cut, as 
we have done in this paper. 

Nevertheless, the very nature of our results suggest a 
fundamental underlying difference between the accretion 
modes we have labeled "cold-mode" and "hot-mode," as 
the two accretion mechanisms show distinct properties 
(even before first infall to the halo virial radius) in terms 
of spatial position and angular momentum content upon 
entering the halo. Whether or not the cold streams break 
up on their way to the disk, the fact that they enter with 
high angular momentum is likely a robust physical dis- 
tinction. Based on the results presented here, we spec- 
ulate that it may in fact be more physically meaningful 
to consider the bi-modal mechanism for gas accretion 
onto galaxies not as fundamentally "cold" versus "hot," 
(though this certainly appears to be true), but rather as 
"filamentary" versus "non-filamentary" modes of gas ac- 
cretion onto galaxy halos. We admit, however, that the 
robustness of any such reclassification would require a 



much more thorough analysis of filamentary versus non- 
filamentary accretion than what we have presented here, 
and would be a useful topic of further study. 
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